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We have investigated the quantum phase transitions in the ground states of sev- 
eral critical systems, including transverse field Ising and XY models as well as XY 
with multiple spin interactions, XXZ and the collective system LMG models, by 
using different quantumness measures, such as entanglement of formation, quantum 
discord, as well as its classical counterpart, measurement-induced disturbance and 
CHSH-Bell function. We have searched for quantum information theory based new 
indicators of the critical points of phase transitions in these critical systems, such 
as measurement-induced disturbance and CHSH function, and have compared and 
contrasted them with the well known critical point detectors, the so-called, quantum 
discord, entanglement and classical correlations. Measurement-induced disturbance 
is found to detect the first and second order phase transitions present in these crit- 
ical systems, while, surprisingly, it is found to fail to signal the infinite-order phase 
transition present in XXZ model. Remarkably, CHSH-Bell function, although is 
not violated, is found to detect all the phase transitions, even when quantum and 
classical correlations are zero for the relevant ground state. 
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I. INTRODUCTION 



A quantum phase transition (QPT) is the qualitative change of the ground state of a 
quantum many body system as the non-thermal parameters (tuning parameters) of the 
system Hamiltonian is varied [T]. QPTs occur at low temperatures, mostly at T = 0, 
where the quantum fluctuations are more dominant than the thermal fluctuations so that 
no classical phase transition is allowed to occur. The superfluid-Mott insulator transition 
in a Bose-Einstein condensate [2] , superconductor-insulator transitions [3] and the quantum 
phase transition in spin models at T = jl] are the remarkable examples of the QPTs. The 
order of the QPT can be determined by looking at the non-analyticities of the derivatives 
of the ground state energy. If the first derivative of the ground state energy with respect 
to the tuning parameter presents a finite discontinuity then the QPT is of the first order. 
Similarly, if the first derivative of the ground state energy is continuous, while its second 
derivative exhibits a finite discontinuity or divergence, then the QPT is second order (also 
called continuous QPT). On the other hand, if any finite-order derivatives of the ground 
state energy is a continuous function with respect to the relevant tuning parameter, then 
the QPT is called an infinite order. 

Quantum correlations present in the ground state of the system are expected to reflect 
the radical change in the ground state of the system at the critical point (CP) of QPTs while 
properly tuning the system's Hamiltonian, because QPTs are related with the change in the 
long-range order of the system. Therefore, one would expect that the signatures of QPTs 
do exist in the measures of quantum correlations. In fact, different measures of quantum 
correlations, such as entanglement [HHS], quantum discord (QD) [TOl - HT] . fidelity [TBI [19] and 
Bell inequalities [201422] have been shown to provide the signatures of QPTs. Although it is 
not a universal role, the extremal points (minimums or maximums) as well as the behavior 
of the derivatives of the correlation measures can signal the existence of a critical point. In 
general, at the first-order QPT point, QD and entanglement (measured by concurrence and 
entanglement of formation (EoF)) are found to be discontinuous at the critical point [TUl 
[131 [H], while at the second-order QPT, the divergence or discontinuity behaviors of the first 
and second derivatives of QD and entanglement could spotlight the critical points [T01T12"] . 
On the other hand, QD and entanglement are found to show local extreme points at the 
CPs of infinite-order QPTs [IMS]. 
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Entanglement may fail to measure the quantum correlations of a state when it reaches 
to a certain level of mixture, especially for distant spin pairs in a spin system and at finite 
temperatures. In addition, QPTs might involve transitions from unentangled to unentangled 
ground states. On the other hand, quantum discord, mostly, does not vanish for mixed 
separable states, actually, almost all quantum states have non-zero quantum discord [23J. 
This feature makes QD more reliable and useful measure compared to entanglement to signal 
the QPTs. It was shown that QD is able to detect the critical points of QPTs, especially 
between far neighbors and also at finite temperatures where entanglement is fully absent 
near the critical points [T2HIH US]. It should be pointed out here that the Bell nonlocality 
as measured by CHSH inequality or its derivative has been used to detect the quantum 
phase transitions only in a few model systems, including Heisenberg XXZ [21] and Kitaev- 
Castelnovo-Chamon models [22]. In these studies, it was shown that CHSH-Bell inequality, 
although is not violated for some cases, exhibits the existence of QPTs present in the models, 
even for separable states. Moreover, in Ref . [20] , the authors investigated only the behavior 
of the CHSH-Bell function near the Ising-like critical point in infinite XY model and found a 
non-ergodic behavior of the nonlocality measure near the phase transition point. However, 
making a general conclusion about the ability of CHSH inequality to detect the CPs of 
QPTs is hard. Although a number of studies indicate a strong relation between the QPT 
and a discontinuity in QD and entanglement or in their derivatives, so far there is no single 
correlation measure that shows universal behavior at the critical points of QPTs. 

The main motivation of the present study is to search for new quantum information 
theory based indicators of quantum phase transitions, such as measurement-induced distur- 
bance (MID) and CHSH-Bell function, and compare them with the well known CP-detectors, 
quantum discord (QD), classical correlations (CC) and entanglement. Towards that goal, 
we will investigate different correlation measures, such as EoF, QD and its classical coun- 
terpart as well as the CHSH-Bell function and the quantumness measure, the so-called, 
measurement-induced disturbance (the explicit definitions of each correlation measure will 
be provided later in the text) in the ground states of several critical systems, including 
transverse field Ising and XY models as well as XY model with multiple spin interactions, 
XXZ model and the collective system LMG model. The focus of the present study is two 
fold. On one hand, we compare and contrast the behaviors of different quantumness mea- 
sures as function of the tuning parameters of the system's Hamiltonians, especially for the 
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regions close to the critical points of the QPTs of the system. On the other hand, we study 
different critical systems that exhibit first and infinite order QPTs as well as several contin- 
uous phase transitions belonging to different universality classes to see whether or not the 
corresponding correlation measures (particularly MID and CHSH-Bell function) can signal 
the abrupt change of the ground state. We have assumed both finite and infinite chains as 
well as the finite size effects of the derivatives of the correlations at the second order QPT 
points. An analysis of QD and EoF for these critical systems has been recently reported 
in Refs. PUJ H21 HE]. Here we both extend and review the results in Refs. [TUl H21 US], and 
also introduce MID and CHSH-Bell function as new theoretical CP-detectors for these sys- 
tems. The major findings of the present study can be summarized as follows: CHSH Bell 
inequality, or its derivative is found to be able to detect all CPs of the QPTs present in the 
considered systems. On the other hand, although MID can signal the first and second order 
QPTs, it fails to detect the infinite order QPT in the XXZ model. Remarkably, in the LMG 
model, the first- derivative of the CHSH-Bell function is found to exhibit the continuous QPT 
between different modes where the other quantum correlation measures as well as CC are 
zero, i.e., they are unable to detect the QPT. 

The paper is organized as follows. In Sec. [TTJ we review the correlation measures for 
bipartite quantum systems and then give the explicit analytic expressions for general X- 



states (a class of states includes Z 2 symmetric states). In Sec. Ill A, we have analyzed the 



quantum criticality in Ising, XY and XY with multiple spin interaction models. In Sees. Ill B 



and III C , we have focused on the ground state correlations and criticality in XXZ and LMG 



models, respectively. A brief summary of important results are given in Sec. IV 



II. CORRELATIONS IN BIPARTITE QUANTUM SYSTEMS 

In this section, we will firstly review the correlation measures for bipartite quantum 
systems and then we will give their analytic expressions in terms of density matrix elements 
for a general two-qubit X-structured density matrix. 

For the model systems that we will consider in the next sections, the Hamiltonians exhibit 
Z2 symmetry; the Hamiltonians are invariant under 7r rotation around a given spin axis 0, 
[TO] . This symmetry leads to the two qubit reduced density matrix at sites i and j in the 
standard basis |1) = |tt) , |2) = III) , |3) = lit) and |4) = \U) (where |t) and \i) are the 
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Pij 



(1) 



eigenstates of the Pauli spin 1/2 ^-operator) has an X-structure [8| [TO] : 

' Mil Uu \ 

U22 U23 
«32 M33 

y U41 M44 j 

For such states, the correlation measures can be calculated analytically. 

We will use entanglement of formation (EoF) as a measure of entanglement which is 
defined as [23] 



EoF( Plj ) = H 



(2) 



where H[x] = —x log 2 x — (1 — x) log 2 (l — x) is the binary entropy function and C(/?y) is the 
concurrence, an entanglement monotone [21]. .EoF quantifies a class of quantum correlations 
that cannot be created by local operations and classical communication only j2UEH]. EoF 
is equal to zero for separable states and one for maximally entangled (Bell) states. For the 
density matrix ([I]), the concurrence can be obtained easily 

C(pij) = 2max{0, \u u \ ~ ^22^33, l^l - y/uuUu}- 

On the other hand, there exist more general quantum correlations that may present even 
for separable states. Recently, Ollivier and Zurek introduced the notion, called quantum 
discord (QD), as a capture of all nonclassical correlations between bipartite quantum sys- 
tems [26]. QD is defined as the difference between the quantum versions of two classically 
equivalent definitions of mutual information and can be given as 



QDfaj) = Ifa) - CC(p 



(3) 



where I(pij) = S(pi) + S(pj) — S(pij) is the quantum mutual information and measures 
the total correlations. Here pij = Trj^pij are the reduced density matrices and S(p) = 
—Tr(plog 2 p) is the von-Neumann entropy. CC(pij) in Eq. ^ is the measure of classical 
correlations between two subsystems [2Z]. It is defined as the maximum information about 
one system that can be obtained by performing a set of projective measurements on the 
other subsystem. Its calculation is based on complex maximization procedure and obtaining 
an analytic expression of CC(pij) for general states is not an easy task. However, for X- 
states the analytic expression of CC(pij) is available if we restrict ourselves to the projective 
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positive valued measurements performed locally on the subsystem j [28] . According to the 
results in Ref. [2E], QD and CC are given as: 

QD( PtJ ) = mm{Q u Q 2 }, CC( Pij ) = ma,x{CC l7 CC 2 }, (4) 

where CCj = H[u u + u 22 ] - Dj and Qj = H[uu + u 33 ] + J2t=i ^ k 1°S2 ^ k + Dj with X k 
being the eigenvalues of pij and H[x] is the binary entropy defined above. Here D\ = H[r] 
where t = (l + - 2(u 33 + u u )] 2 + 4(|w 14 | + |m 23 |) 2 ) /2 and D 2 = - Yl=i u kk log 2 u kk - 
H[u\i + U33]. In general, QD or CC are not symmetric quantities, i.e., their values can 
depend on which the projective measurements are performed on. Indeed, for general X 
states, the order of measurements performed are important since S(pi) 7^ S(pj). However 
for the considered model systems in the present study (except LMG model), m 22 = u 33 which 
ensures S(pi) = S(pj); irrespective of the measurement is performed locally on subsystem i 
or j. 

Along similar lines, the so-called measurement-induced disturbance (MID) was proposed 
very recently by Luo in order to quantify the quantumness of correlations [29]. The idea on 
the definition of MID is based on that for a given state p, any complete local projective 
measurements {n" £8>n^} on the subsystems a and b will induce a classical state n(p). If the 
measurements {n^} are chosen as the spectral resolutions of the reduced states; p a = J2 i p^n° 
and pb = J2jP b i^ b j (where p^ ,b are the probability amplitudes), the state U(p) = J2ij(J^i ® 
n^)p(n" cg> n^) will be the closest classical state to p. Indeed, the measurements induced by 
the spectral resolution of the reduced states will leave the marginal states invariant; in the 
sense of the least disturbing measurements. MID is just the quantum mutual information 
difference between the given state and the closest classical state to p: 

MID = I(p)-I(Tl(p)). (5) 

MID is easily computable and can be applied to bipartite higher dimensions. For the density 
matrix in Eq. (fXT) , the closest classical state would be U-(pij) = diag{un,u 2 2, U33, u^}- Since 
the reduced states are not affected by this measurements and the closest classical state is 
diagonal, MID reduces to the entropy distance, S(x\\y) = —tr(x\og 2 y) — S(x), between 
the given state and the closest classical state to pif 

MID(p tJ ) = S( PlJ \\U(p v )) 

4 4 
= ^2\ k log 2 X k - ^ u kk log 2 u kk , (6) 

k=l k=l 
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where \ k are the eigenvalues of py. 

The last correlation measure that we will use in the present study is the CHSH-Bell 
inequality [30]. The violation of CHSH inequality describes the part of quantum corre- 
lations (often called nonlocal correlations) that cannot be reproduced by a local hidden 
variable model, i.e. by classical systems. Although entangled pure states violate the CHSH 
inequality, mixed entangled states may not imply that their correlations cannot be classi- 
cally reproduced with certainty. By using the Horodecki criteria [31], the maximum of the 
CHSH Bell function for bipartite X-states can be obtained easily [32] : 

CHSH{p i3 ) = max{ J B 1 , B 2 }, (7) 

where B\ = 2^/k\ + Ak\ and B 2 = 4y7c| + fcf, k x = (mh+M44-«22-«33), k 2 = (\uu\ + \v^3\) 
and k 3 = (\uu\ — | w 2 3 1 ) • When the value of CHSH function exceeds 2 for a given state , it 
is violated, that is no classical local model can reproduce all correlations of these states. 

The correlation measures (EoF, MID and QD) are equal to each other for bipartite pure 
states, but they disagree on the value of correlations for mixed states. MID and QD may 
give non-zero values for separable states. On the other hand, MID can be nonzero, even 
maximal for purely classical states (where QD = 0). Therefore, MID is sometimes regarded 
as an unfaithful and nonrefined measure of nonclassical correlations in bipartite states [33- 
135] . CHSH Bell function attains its maximum value for maximally entangled states, but in 
general it is violated for states having high entanglement. 

III. MODEL SYSTEMS, QUANTUM PHASE TRANSITIONS AND GROUND 

STATE CORRELATIONS 

A. XYT model and its variants 

The anisotropic spin 1/2 XY chain with three spin interaction under external transverse 
field (XYT model) can be defined by the following Hamiltonian [16, 36] : 

hxyt = -t{ (^) + i^f) ^ + x ° z ] 

N 

»K • (8) 

i=i 
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where a* ,v ' z are the spin 1/2 Pauli matrices on site j, N is the total number of spins in the 
chain, 7 is the degree of anisotropy and A and a are the strengths of the external field and 
the three spin interaction, respectively. 

The XYT model reduces to the transverse field Ising model for a = and 7 = 1 and 
to the anisotropic XY model for a = 0. It is exactly solvable and can be diagonalized by 
mapping it into spinless free fermion model with single orbitals. This can be implemented 
through Jordan Wigner map followed by a Bogoliubov transformation (see, e.g., Refs. [T61I36] 
for the diagonalization of XYT Hamiltonian). 

The symmetries of the Hamiltonian of the system simplifies the structure of the two point 
reduced density matrix; Since Hamiltonian of XYT model is real, then pij = p*j, because 
of translational invariance of Hxyt, the two-point density matrix of the system is position 
independent, i.e., p^ = pu+ r for integer r. Furthermore, the density matrix commutes with 
of crj because of the global rotation symmetry (^-symmetry) of Hxyt- As a result of these 
symmetries, the two qubit reduced density matrix of the system has X-structure with the 
following matrix elements [16J: 

«n = \ (1 + 2 (a*) + (ofoU,)) , 
U22 = u 33 = - (l - (ofof +r » , 

U23 = U 32 = - «« r ) + (Oi<>i+r)) , 
Uu = U41 = 1 «0f 0f +r ) - (of of+r)) > 

uu = \ (1 - 2 (cr z ) + (aM +r )) . (9) 

In the following, we will restrict ourselves to the ground state correlations (where the 
temperature T = 0) between the nearest-neighbor sites (r = 1). The transverse magnetiza- 
tion (a z ) and the two-point correlation functions between nearest-neighbor sites at T = 
for the XYT model can be calculated as [37, 38J: 

(ofof ¥1 ) = G- 1 , <afaf +1 )=G 1 , 

(alal +1 ) = {a z f - G X G- Xl (10) 
where G k = - J2k (cos(x k r)e k + 7 sin(x fc r) sin(x k )) /Nu k and uj k = ^ e\ + (7 sin(x fc )) 2 , e k = 
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A — cos(xfc) — 2acos(2xfc) is the energy spectrum with Xk = ^r, k = —M,...,M and 
M = (N - l)/2 for odd N. 

Now, we discuss the critical behavior of correlation measures for the transverse field Ising, 
XY and XYT models. 

I. Transverse field Ising model 

As a first investigation, we set a = and 7 = 1 in the XYT Hamiltonian which reduces 
to the well known transverse field Ising Hamiltonian [To] : 

N 

II, ,,,,, »;<., • Arrj}. (11) 

i=i 

The system will undergo a second-order quantum phase transition at the critical value A c = 1 
where the global phase flip symmetry breaks and the correlation length diverges [3J. Our aim 
in this section is two-fold; first we want to analyze the behavior of correlation measures (EoF, 
QD, CC, MID and CHS H-Bell function) near the critical point. Secondly, we want to 
check the ability of the correlation measures to detect the QPT and investigate the finite-size 
scaling behaviors of the derivatives of the correlations at the critical point, A c . 




0.0 0.5 1.0 1.5 2.0 

A 

FIG. 1. (Color online) Correlation measures as a function of A for nearest-neighbor spins in the 
transverse field Ising model (7 = 1 and a = 0) for a finite chain with N = 2001 sites. 

In Fig. 1, we display the correlation measures as a function of A for a chain with N = 2001 
sites in the transverse field Ising model. When the parameter A = 0, all spins are oriented 
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Pii+1 



(13) 



in either positive or negative z-direction and the ground state is doubly degenerate given by 
the density matrix: 

P= \ 1+ + •••+) + + 1 - |- -•••-) (---..-I, (12) 

where |±) = (|t) ± \\))/ \/2. After taking a partial trace of p overall spins except the nearest 
neighbors, the two-point reduced density matrix in the standard basis would be 

/ 0.25 0.25 \ 

0.25 0.25 

0.25 0.25 

\0.25 0.25 / 

which contains EoF = QD = 0, but MID = CC = 1 as can be seen from Fig. 1. On the 
other hand, for A — > oo all spins are oriented along the positive z-direction and the ground 
state density matrix is equal to p = |tT • • • T) (tt • • • Tl which is a product state that contains 
no type of correlations between any parties. The transition between the paramagnetic and 
ordered ferromagnetic phases occurs at A c = 1. 

As the magnitude of external field increases from zero, QD and EoF increase and they 
attain their maximum value near the critical point [TQl[TS]. On the other hand, CC and MID 
are monotonous functions of A; the magnitudes of CC and MID decrease as A increases. 
Similar to EoF and QD, CC and MID increases when approaching the critical point from 
right for A > 1. CHSH Bell function, although below 2, reaches its minimum value in a 
region that is very close to A c = 1 (see inset in Fig. 1). Fig. 1 also shows that A-dependence 
of MID, peculiarly, is very close to that of CC and both of them are equal to 1 for A = 
where EoF = QD = 0. One should note that although MID is easier to calculate than QD, 
it shows an inconsistency to reflect the quantum correlations compared to QD. Its being a 
true quantum correlation measure as defined in Sec. II, is questioned, recently [33H35] . It 
is claimed that lack of any optimization procedure over the projective measurements in the 
definition of MID makes it a weak correlation measure for some states [3"31 - f3"5] . 

The QPT at A c = 1 is a second-order phase transition, thus only the derivatives of the 
correlation measures can signal the presence of the QPT as can be seen in Fig. 2. Indeed, the 
first derivative of EoF (Fig. 2(a)), CC (Fig. 2(c)), MID (Fig. 2(d)) and CHSH (Fig. 2(d)) 
show pronounced minimums or maximums at some values of A which tends to the quantum 
critical point, A c = 1, as the system size approaches infinity (N — > oo). The insets in Fig. 2 
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FIG. 2. (Color online) The first derivative of correlation measures ((a)-(e)) and the second deriva- 
tive of QD (f ) for nearest-neighbor spins as a function of A in the transverse filed Ising chain (7 = 1 
and a = 0) for lattice sizes N = 51 (blue-dotted), N = 101 (green-dotdashed), ./V = 501 (red- 
dashed) and N = 2001 (black-solid). Insets: Finite size scaling behaviors of the derivative(s) of 
the correlations at the critical point A c = 1. Here filled circles are the numerical data, while the 
lines are the fitting curves fitted to y = ax + b for (a)-(e) and to y = ax 2 + bx + c for (f). 



show that these correlations exhibit a logarithmic divergence at A c . Peculiarly, the first 
derivative of QD behaves differently; it shows an inflection point around A c (see Fig. 2(b)), 
while its second derivative shows a pronounced maximum which has a quadratic logarithmic 
divergence at A c (see Fig. 2(f) and Ref. [10]). It is remarkable to note that the behavior of 
QD is rather different from the other correlation measures considered in this paper. 
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2. Transverse field XY model 

The second model that we will consider in this study is the transverse field XY model 
which can be obtained by setting a = in Eq. (JsJ) : 

Hxy = -J2{ (^) -H+i + (V) ^ + ■ (14) 

The XY model has an interesting quantum phase transition belonging to two different 
universality classes [39J: (i) For non-zero values of anisotropy 7, the model exhibits a second- 
order quantum phase transition at the critical point A c = 1. In the interval < I7I < 1, 
this quantum phase transition has a universality class which is the same as that of the 
transverse field Ising model, (ii) For A G (0,1), the model exhibits again a second order 
phase transition, the so-called anisotropy transition, at the critical point 7 C = 0. This critical 
point separates a ferromagnetic ordered along x-direction and the ferromagnet ordered along 
the ^-direction. 

In Fig. 3, we display EoF (a), QD (b), CC (c), MID (d) and CHSH -function (e) for 
nearest-neighbor spins in the XY spin chain with N = 2001 as a function of 7 and A at 
zero temperature. All figures demonstrates that there is a pronounced difference in the 
correlation measures between the regions < A < 1 and 1 < A < 2 where the QPT occurs 
at the critical line A c = 1. When A < 1, EoF and QD decreases with 7. Contrarily, CC 
and MID increases. On the other hand, when A > 1, the correlations (MID, QD, EoF 
and CC) increase with anisotropy. Actually, MID and CC are increasing functions of 7 
independent of A. C H S H-Be\l function is always below 2, i.e., Bell inequality is not violated, 
and its dependence on 7 is weak for the regions far from the critical point A c = 1. Near 
the critical point, CHSH attains its minimum. The minimum decreases when anisotropy 
decreases and becomes more sharp when 7 — > (XX model). 

As mentioned above, the XY model undergoes second-order Ising like quantum phase 
transition at A c = 1. Thus, in Fig. 4, we display the first derivatives of correlations as a 
function of A for different values of 7 and N = 2001. As depicted in these sub figures, the 
first derivatives of all correlation measures are divergent at A c = 1, signaling the presence 
of the QPT. It seems that the decrease in 7 sharpens the singularities to some degree. One 
should note that the derivative of QD has a different character for XY model compared to 
that for the Ising model discussed above. For finite chain Ising model (7 = 1 and N = 2001), 
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(b) QD 




(d) MID 

0.0 




FIG. 3. (Color online) EoF (a), QD (b), CC (c), MID (d), CHSH (e) in the XY spin chain (a = 
0) with N = 2001 sites as a function of 7 and A. 

the first derivative of QD is found to exhibit an inflection point around A c , while its second 
derivative is divergent at A c (see Figs. 2(b) and 2(f)). As shown in Fig. 4(b), for anisotropic 
XY model, its first derivative is already divergent at A c . 

In Fig. 5, we plot EoF, QD, MID, CC and CHSH-Bel\ function (inset) as a function of 
anisotropy 7 for iV = 2001 and A = 0.1. All correlation measures are symmetric according 
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FIG. 4. (Color online) The first A-derivative of correlations, EoF (a), QD (b), CC (c), MID (d) 
and CHSH (e) as a function of A in the anisotropic XY model (a = 0) for different values of 
anisotropies 7 = 0.3 (black-solid), 7 = 0.7 (blue-dotted) and 7 = 0.9 (red-dashed) and N = 2001. 



to 7 —7. EoF and QD are maximum at 7 C , while MID, CC and CHSH reaches their 
minimums. It is safe to say that all the considered correlation measures precisely detect 
the anisotropy phase transition by reaching their maximums or minimums. Especially, 
the behaviors of QD and CC are more drastic, showing extremum at 7 C with a cusplike 
behavior arising from the minimization procedure over the complete set of projectors in 
their definitions. In fact, the first 7 derivative of EoF, MID and CHSH show an inflection 
around the critical point 7 C = 0, while the first derivatives of QD and CC are discontinuous 
at 7 C . On the other hand, second derivatives of all correlations diverge at the anisotropy 
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y 

FIG. 5. (Color online) QD (red-dashed), EoF (blue-dotted), CC (green-dotdashed) , MID (black- 
solid) and CHS ff-Bell function (inset) as a function of 7 for the anisotropic XY model (a = 0) 
with N = 2001 and A = 0.1. 

transition point (not shown here). 

3. XYT model 

Now, we consider the general Hamiltonian, Eq. (JsJ) , and analyze the correlation measures 
in the ground state of the XYT model. The model has a rich quantum phase diagram (see 
Refs. [521 ED])- Differently from the transverse field XY model, the QPTs are caused not 
only by the external field and anisotropy, but also by the three spin interactions. Indeed, 
the competition between A, 7 and a may induce QPTs belonging to different universality 
classes, for example, ferromagnetic, chiral and different spin liquid phases [36| HO]. 

We start the discussion of the correlation-phase diagram of XYT model with an investi- 
gation of the effect of three-spin interaction strength and the external magnetic field while 
the anisotropy is held constant. In Fig. |6j we display EoF and MID as a function of A and 
a at 7 = 0.5 and N = 2001 in the XYT model. The dramatic changes in EoF and MID 
show the phase transitions in the model (similar for other considered correlation measures, 
thus we have not displayed them here). The region between a G [0,0.5) and A G (0.9,2.0) 
indicates the spin-saturated phase (ferromagnetic), while the regions between a G (0.5, 1.5) 
and A G [0,2] as well as a G (1.5,2] and A G [0,2] are the two kinds of incommensurate 
phases with weak spontaneous magnetization (see Refs. [TBI EH HI] for details). 

All transitions are second-order, thus one can expect that the derivatives of correlation 
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FIG. 6. EoF (a) and MID (b) as a function of A and a at 7 = 0.5 and N = 2001 in the XYT 
model. 

measures could signal these QPTs. In Fig. [7| the a-derivative of correlation measures as a 
function of a are displayed at A = 2.0 and 7 = 0.5 in XYT model with different lattice sizes. 
The first derivatives of all correlation measures diverge at a c ,i = 0.5 and a Cj 2 = 1.5 showing 
the presence of QPTs. These extremum points indicate that there exits three different phases 
as function of a in the region a G [0, 2] for the XYT model. Note that the first derivatives 
of all the correlations exhibit logarithmic divergence at a c l and a Cj2 as N approaches to 
infinity (see insets in Fig. [7]). 

As shown in Fig. [6j for a fix a, the external magnetic field can also lead to QPTs in 
XYT model. For instance, for a = 0.2 and a = 1.0, there are second-order quantum 
phase transitions at A c = 1.4 and A c = 1.0, respectively. The divergent behavior of the 
first derivative of correlations with respect to A signals these phase transitions as shown in 
Figs. [8]^a) and ^b) . They diverge logarithmically as the system size increases, which is 
similar to the QPTs in the line of a (we have not displayed them here). In Ref. [TBj, the 
authors studied QPT in XYT model and they only focused on the critical point a c = 0.5 
for 7 = 0.5 and A = 0. It was shown that the first a-derivative of QD and CC are capable 
detecting the QPT. Different from Ref. |16j . we have shown that QD and CC as well as 
MID, CHSH and EoF can signal all the phase transitions present in the XYT model for 
7 = 0.5 and A > 0. 
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FIG. 7. (Color online) The a-derivative of correlation measures as a function of a at A = 2 and 
7 = 0.5 in the XYT chain with iV = 51 (blue-dotted), N = 101 (red-dashed) and N = 2001 (black). 
Insets show the finite size scaling behaviors of correlation measures at the QPTs a c ,i = 0.5 and 
a C) 2 = 1.5 where the filled circles are numerical data, while the lines are the fitting curves fitted to 
y = ax + b. 

B. XXZ model 



The next model that we will tackle in the present paper is the one dimensional anisotropic 
Heisenberg XXZ chain given by the Hamiltonian 



N 



H 



XXZ 



E 5 ?s 



f ' +S^Sf +1 + AS!Sf 



i+l 



(15) 



i=i 
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FIG. 8. (Color online) First A-derivative of correlation measures as a function of A in the XYT 
model with N = 2001, 7 = 0.5 and a = 0.2 (a) and a = 1.0 (b). 



where A is the anisotropy parameter, Sj = <jf/2 and N is the number of particles in 
the chain. The model has two critical points that separates three different phases: (i) At 
A c = —1, there is a first order QPT that separates the ferromagnetic phase, where all the 
spins are pointing in the same direction for A < — 1, and the gapless phase, where the 
correlation decays polynomially for —1 < A < 1. (ii) At A c = 1, there is an infinite order 
QPT that separates the gapless phase and the anti-ferromagnetic phase for A > 1. 

The Hamiltonian (15) has spin-flip symmetry that leads to zero magnetization (a z ) = 
and it is symmetric with respect to to the rotation in xy plane which leads to (erferj) = 
(afa^j. Due to these symmetries, the nearest neighbor reduced density matrix would be in 
X-state form as in Eq. ([T| with matrix elements 



Mil = «44 
U23 = U32 



4 
4 

<«fi) 



Ul4 = 0. 



(16) 



The XX Z model cannot be diagonalized, but the two point correlation functions at zero 
temperature, T = and in the thermodynamic limit, iV — » 00 can be obtained by using 
Bethe ansatz technique [42]. The spin-spin correlation functions between nearest-neighbor 
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spin sites for — 1 < A < 1 are given by [13] 

2 [°° dx x cosh x 2cot(7r$) 



7r 2 J_ OQ sinhx cosh 2 ($x) 
cos7r<3? f°° dx x cosh x 



7T 
1 



dx sinh ((1 — $)a;) 
sinh x cosh $s 
fix sinh ((1 — 



7r 2 y.oo sinh x cosh 2 vr sin(7r<£>) J_ OQ sinhx cosh$x 
where = - cos -1 A. For A > 1 the correlation functions are given by [H] 

"°° +i/2 dx 



(17) 



(aM + i) = l + 2 



oo+i/2 sinh(7rs 

°° +i/2 rfx / r 



cot(vx) coth(z/ 
- cosh v 



x 



sm (z/a;^ 
cot(^x) 



-oo+i/2 sinh(7rx) \sin 2 (z/x) ~ sinhz/ 
where v = cosh -1 A. And for A < —1, (of = 1 and (crfa-f +1 ) = 



(18) 




FIG. 9. (Color online) (a) shows EoF (red-dashed), QD (black-solid) and CHSH (inset), while 
(b) shows CC (red-dashed) and MID (black solid) for nearest-neighbor spins as a function of A 
in the XXZ chain. 

In Fig. [9j we display the correlation measures as a function of anisotropy parameter A 
in the XXZ chain at T = and N — > oo. At the first order critical point, A c = —1, all 
correlation measures, except EoF, are discontinuous showing the presence of QPT, while 
this jump behavior is not observed in EoF. At the infinite order QPT, A c = 1, QD, CC 
and CHSH function show abrupt cusps. The cusp behavior in the considered correlation 
measures is due to the level-crossing in the model which involve transition from the gapless 
to antiferromagnetic phases. On the other hand, EoF is less sensitive to this level crossing 
it just reaches its maximum at A c = 1 which is also the case for QD [HJ [15]. These 
cusplike behaviors at infinite order QPT point implies that the first derivatives of QD, 



20 



CC and CHSH are discontinuous at A c = 1, while their second derivatives are divergent. 
Surprisingly, MID (black line in Fig. (9^b)) does not exhibit any abrupt cusp or extremum 
point at A c = 1. It is continuous at the phase transition point and shows no sign of the 
infinite order phase transition. By using the correlation functions in Refs. [H2 E]> we have 
also calculated MID for next-nearest neighbors, and have observed that MID cannot signal 
the infinite order phase transition for this case as well. It is remarkable to note that between 
< A < 1, QD and EoF increase with A, while MID is decreasing function of anisotropy 
in this region. On the other hand, they are all decreasing functions of A for A > 1 and 
reaches zero for A — > oo. We should note that while MID and CC behave similarly for 
Ising and XY models as functions of tuning parameters, as discussed previously, they exhibit 
completely different dependence on A for the XX Z model. See also Ref. [21 J, where the 
authors showed that CHSH-Bell function, although is not violated, is found to capable to 
detect the presence of QPT at A c = 1 for next- and next-next-nearest neighbor spins where 
entanglement is absent for such distances in the XX Z model around the critical point. 



C. LMG model 



Finally, we will analyze the ground state correlations of the model introduced by Lipkin, 
Meshkov and Glick {LMG) jl5]. The LMG model, originally, describes a two-level Fermi 
system with one-level below the Fermi level, while the other level just above. The level 
below the Fermi level is assumed to be filled with V nucleons (here V can also be thought 
of as degeneracy). In its spin 1/2 representation, the LMG model can be regarded as a 
chain of N spin particles with infinite-ranged interactions, where each spin is subject to an 
external transverse magnetic field A. Indeed, the Hamiltonian of the LMG model can be 
written as [9j M, HE] : 

Hlmg = — {Sx — Sy) + AS^, (19) 

where S z = X/m=i I ( c +m c +m ~ cL m c_ m J and S x + iS y = Y2 m =i c +m c -m, with c^ m and c^_ m 
being the particle creation operators in the upper (+m) and lower (— m) levels, respectively. 
The model has a second-order quantum phase transition at A c = 1 between a symmetric 
phase (A > 1) and a broken phase (A < 1) [10J H?]. The ground state is unique and fully 
polarized in the direction of the field for A > 1, while its is twofold degenerate for A < 1. 
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The ground state of the LMG model for N — > oo can be obtained by using Hartree- 
Fock (HF) approach. Indeed, HF ground state provides the exact description of the critical 
point. The pairwise ground state density matrix for general modes (+m) and (— n) in the 
HF ground state of the LMG model can be found as pH] 



/ 



Pm,—n 











(M +m M_ n ) 

{M +m N_ 



\ 















{N +m M_ n } 

(iV+ m AL n ) J 



(20) 



where = 1 — and = c k Ck and the matrix elements are given as 



(M +m M_ n ) 



)sin 2 a cos 2 a, 



(M +m N_ n ) = 5 mn cos 2 a + (1 - 5 mn ) cos 4 a, 



(N +m M_ n ) 
(N +m N_ n ) 



C+mP—n 



5 mn sin a + (1 - 5 mn ) sin a, 



'1 — 5 mn )sin 2 a cos 2 a, 



5 mr) cos a sin a, 



(21) 



where 5 m „ is the Kronecker-Delta function and a 



| cos 1 A for A < 1 and a = for 



A > 1 (see Refs. pU] and [46J , for the detailed calculations of Eqs. (20) and (21)). One 
should note that the exact ground state of the LMG model, instead of the HF ground state 
studied here, can be obtained and highly entangled for A > 1 [38] . 



In Fig. 10, we display the calculated correlations and their first derivatives with respect to 
A between the modes (+m) and (— m) as well as between the modes (+m) and (— n) (m ^ n) 
in the HF ground state of the LMG model. CC, QD, MID and EoF are all equal for 
the density matrix of the ground state for the (+m, —m) modes, because for those modes 
the ground state is pure and these measures are equal for the pure states. The value of 
correlations range from 1 for A = to zero for A > 1 (see Fig. [To^a)). Also CHSH Bell 
function is violated for A < 1 (Fig. [To](b)). The derivatives of correlation measures are 
discontinuous at A c = 1, signaling the QPT (see insets in Figs. [To|(a) and [lO](b)) . On 
the other hand, for m ^ n, there is no correlations (MID, QD, CC and EoF) between 



the modes. In fact, the density matrix Eq. (20) is diagonal for m ^ n and contains no 
pairwise classical and quantum correlations. However, the diagonal elements vary with A 
and the pairwise density matrix changes between maximally mixed state for A = and the 
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FIG. 10. Correlations (QD = EoF = CC = MID) (a) and CHSH-Bell function (b) between 
the modes (+m) and (— m) as a function of A in the HF ground state of the LMG model, (c) 
shows CHSH function for the modes m ^ n. Insets are the corresponding first order A-derivative 
of correlation measures as a function of A. Note that for m 7^ n, there is no bipartite correlations, 
i.e., QD = EoF = CC = MID = 0. 



fully polarized state for A > 1. Although this transition does not contain any bipartite 
correlations, remarkably, CHSH Bell function changes with A and the discontinuity (jump) 



in its first derivative at A c = 1 exhibits the QPT (see Fig. 10 c)). 



IV. CONCLUSIONS 



We have investigated the quantum phase transitions of the critical systems, such as 
transverse field Ising, XY, XYT, XXZ and LMG models by using different quantumness 
measures, such as entanglement of formation, quantum discord, as well as its classical coun- 
terpart, measurement induced disturbance and CHSH-Bell function. Different aspects of 
quantum correlations of a given state are compared and contrasted when the tuning pa- 
rameters of the system's Hamiltonians are changed, specifically at different phase regions. 
The ability of correlation measures to detect different order QPTs as well as the finite-size 
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scaling behaviors of them at different continuous phase transition points have been also 
investigated. 

The most important findings of the present work can be summarized as follows. The 
correlation measures that we consider reflect somewhat different aspects of the quantum as 
well as classical properties of the ground state of the examined models and it is not possible to 
find a universal behavior of these measures for signaling the quantum phase transition critical 
point. Furthermore, for individual measures, the dependence on the tuning parameter near 
the critical points is not universal across the order of the transition. In particular, for the 
first order QPT observed in XXZ model, all the considered measures, except EoF, show 
a discontinuity behavior at the CP, while at the infinite-order QPT observed in the same 
model, MID has no indication of the QPT. On the other hand, EoF has its maximum and 
QD, CC and CHSH have cusp-like structure, showing the infinite-order phase transition. 
The number of second-order QPTs investigated in the present work is much higher than that 
of the first or infinite-order QPTs. None of the considered measures was found to behave 
exactly the same at CP for the investigated seven second-order phase transitions. Except 
for the LMG model for (+m, —n) modes with m 7^ n, the first or second derivatives of 
all correlation measures are found to be discontinuous or divergent at the critical point of 
the relevant models, signaling the CPs of QPTs. Remarkably, the CHSH function is the 
sole indicator of the phase transition between (+m) and (—n) modes for m 7^ n in the HF 
ground state of LMG model. 

The main motivation of this study was to search for new quantum information theory 
based indicators of phase transitions, such as MID and CHSH-Bell function, and compare 
them with the well known CP-detectors, QD, CC and EoF. Remarkably, the CHSH- 
Bell function is found to exhibit all the phase transitions present in the considered critical 
models, even when there is no bipartite classical or quantum correlations. However, one 
should note that the bipartite states (except LMG model for the (+m, —m) mode case) 
are mixed and do not violate the CHSH inequality, although they contain entanglement. 
This indicates that the correlations possessed by the ground states are highly local. On the 
other had, CHSH function is a non-ergodic quantity, i.e., it does not attain its equilibrium 
value, until the state is frozen. We believe that this non-ergodicity behavior of the CHSH 
function makes it more reliable and useful correlation measure to indicate QPTs than the 
other correlation measures considered in this paper. However, its ability to detect the CPs 
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of QPTs in different critical systems and particularly at finite temperatures needs strict 
investigations to make a direct connection between CHSH function as a CP-detector and 
experiments. 

It would be also very interesting to study different quantumness measures (see the recent 
review article [IS]), especially the ameliorated MID to see whether or not this amplified 
version of MID can signal the infinite order phase transition in XXZ model. Also, the 
analysis to see whether or not this is a generic feature of MID at infinite order phase 
transitions exhibited by different models [7J [T7j needs further considerations. 
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